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Abstract — Using the Archeops and WMAP data we perform 
a study of the anti-center Galactic diffuse emissions - thermal 
dust, synchrotron, free-free and anomalous emission - at degree 
scales. The high frequency data are used to infer the thermal 
dust electromagnetic spectrum and spatial distribution allowing 
us to precisely subtract this component at lower frequencies. 
After subtraction of the thermal dust component a mixture of 
standard synchrotron and free-free emissions does not account 
for the residuals at these low frequencies. Including the all-sky 
408 MHz Haslam data we find evidences for anomalous emission 
with a spectral index of -2.5 in Trj units. However, we are not 
able to conclude regarding the nature of this anomalous emission 
in this region. For the purpose, data between 408 MHz and 20 
GHz covering the same sky region are needed. 

Index Terms — ISM: general - ISM: clouds - Methods: data 
analysis - Cosmology: observations - Submillimeter - Catalogs 

I. Introduction 

The anomalous microwave emission (AME in the following), 
is an important contributor of the Galactic diffuse emissions 
in the range fro m 20 to 6 GHz. I t was first ob served 
by (de Oliveira- C osta et all 119971: iKogut et all 119961) 
and then identified by "dLeitch et all I1997D as free-free 
emission from electrons with temperature, Te > lO^K. 
iDraine & LazarianI (Il998al) argued that AME may result from 
electric dipole radiation due to small rotating grains, the 
so-called spinning d u st. Mo dels of the spinning dust emission 
dPraine & LazarianI Il998bl) show an electromagnetic spectra 
peaking at a round 20-5 GHz being able to reprodu ce the 
observations (Finkbeiner, 2003'; "de Oliveira- Costa et alj. 
IWatson et al., 2005; Idesias-Groth, 2005; Casassus et alj. 
12008; Dickinson et al., 2009; Tibbs et al., I2010D . The initial 
spinning dust model has been refined reg arding the shape and 
rotational properties of the dust g rains (|Ali-Haimoud e t al., 2009; 
iHoang et aU .2010. .201k .Silsbee et all 120111). An alternative 
explan ations of AME was proposed bv lDraine~& Lazaria^ 
based on magnetic dipole radiation arising from hot 
ferromagnetic grains. This kind of models associated to single- 
domain predict polarization fraction much bigger than the 
electric dipole ones ( Lazarian & Draine, 2000). Origi nal models 
have been mainly ruled out by many studies (Battistelli et al., 
2006: Casassus et aL, 2006; Kogut, 2007; Mason et al., 2009; 
Lopez- Caraballo et all l201lh although modern variants of those 



may still be of interest (B. Draine private communication). 

Correlation between microwave and infrare d maps, 
mainly dominated by dust thermal emission iDesert et all 
1990), was observed for vario us exper ime nts, for 
exampl e on COBE/ DMR ( Kogut et all Il996dlbh 
OVRO dLeitch et allll997D . Saskatoon (de Oliveira- Costa et al 
19970 . survey at 19GHz ide Oliveira- Costa et al 



Tenerife 
signal was find in small 
and in some molecular 



dde Oliveira-Costa et all 119991). S imilar 



COSMOSOMAS 
2005), AMI 
CBI ( Casassus et al. 



region 

clouds 

/Genova-Santos et al., 
( Ami-Consortium: Scaife et al. 



by dFinkbeineil 120031) 
based ^on data from 

2oTi_ 



Watson et all 



-rT««w l^OOd; iCastellanos et al., 2011), 

VSA dTibbs et all 120101) and Planck ( Planck- Collaborationj. 
2011). Recent studies based on several sets of data dBot et all 



20101) found similar results. 



Independently, ^Be nnett et alJ d2003l) proposed an alternative 
explanation of AME based on flat-spectrum synchrotron 
emission associated to star-forming regions to explain part 
of the WMAP first-year observations. This hypothes is seems 
to be disagreement with results fr om |de O iveira- Costa et al, 
(2004); Fernandez-Cerezo et alJ d200^ fffildebrandt et al 



(2007); Ysa rd et alJ (l2010h which showed that spinning dust 
was the most trusta ble emission t o explain the excess below 20 
GHz. Furthermore, Davi es et alj (12006) showed the existence 
of important correlation between microwave and infrared 
emission in r egion s outside star-forming areas. More recently, 
iKogut et"al] d2011l) discussed the fact that spinning dust fits 
better to ARCADE data (3.8 and 10 GHz) than a fiat-spectrum 
synchrotron. 

We propose here to study the Galactic diffuse emissions in the 
Galactic plane, particularly focusing on the anti-center region. 
The observational data, from 408 MHz to 3000 GHz, used for 
this study are presented in Section [III Section IIIII discusses in 
details the contribution of the diffuse Galactic thermal dust 
emission using the high frequency data. In Section Jvl we consider 
a simple free-free and canonical synchrotron emission model 
for the thermal dust subtracted microwave data. The possible 
contribution from anomalous emission is discussed in Section |Vl 
We draw conclusions in Section IVll 

II. Microwave and millimeter observations 

We describe in this section the data used for the analysis 
presented in this paper. As we are interested in the Galactic 
diffuse emission we consider only large coverage sky surveys in 
the radio, microwave, millimeter and infrared domain including 
the 408 MHz all-sky survey and the Wmap Archeops and Iras 
data. 

408 MHz all-sky survey 

In the radio domain, the 408MHz all-sky continuum survey 
(Haslam et al. (1982)) at a resolution of of 0.85 degrees, is a 
good tracer of the synchrotron emission. In particular, we use 
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TABLE I 

RMS OF THE HIGH FREQUENCY DATA AND OF THE RESIDUALS AFTER SUBTRACTION OF THE DUST MODEL. 



Frequency (GHz) 


Data nns (mK/^j) 


Residual rms (mK/^j) 


Noise standard deviation (mKjij) 


143 


0.194407 


0.0240274 


0.0241575 


217 


0.315175 


0.0306455 


0.0406982 


353 


0.498923 


0.0705655 


0.0404479 


545 


0.699145 


0.113738 


0.157994 


3000 


0.0973817 


0.00830546 


0.0197808 


5000 


0.0134645 


0.00974734 


0.00176314 
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Fig. 1. From top to bottom and from left to right: maps of the best-fit thermal dust emission spectral indices and temperature and uncertainties at 2 a (95 
% C.L.) (right). 



the 408MHz all-sky map available on the LAMBDA website in 
the HEALpix pixelisation scheme (Gorski et al., 2005). The 408 
MHz all-sky survey map was smoothed down to a resolution 
of 1 degree and do wngraded to Nside — 64 in the HEALpix 
pixelisation scheme dGorski et alj d200^ ). The u ncertainties on 
this m ap are assumed to be of 10 % following iHaslam et alj 
dl982h and are mainly due to calibration errors. 

Wmap 

In order to estimate the diffuse Galactic emission at microwave 
frequencies we used the maps in temperature using the K, Ka, Q, 
V and W band maps of the Wmap mission of its 7-years Wmap 
citephinshaw09. In particular, we used the co-added maps 
available on the the Lambda web site, also smoothed down 
to a resolution of 1 degree and downgraded to Nside — 64. 
Uncertainties in the Wmap data were compu ted assuming a 
uncorrelated anisotropic noise as described in (Hinsha w et all 
12009). The variance per pixel at the working resolution was 
computed using the variance of a single hit and the number 
of hits per pixels. 

Archeops 

In the millimet er wavelength s we use the Archeops bal- 
loon experiment dBenoit et all 120041) data. Archeops ob- 



served the sky at four frequency bands: 143, 217, 353 and 
545 GHz, with a resolution of 11. 13, 12 and 18 arcmin 
respectivelv dMacias-Perez et all |2007|) . The Archeops survey 
covers about 30 % of the sky mainly centered in the Galactic 
anti-center region. We use here the original Archeops maps 
which were also smoothed down to a resolution of one degree 
and downgraded to Nside = 64. 



Iras 

In the infrared, we have used the new generation of the Iras 
data (InfraRed Astronomical Satellite) at 100 and 60 /xm (3000 and 
5000 GHz). This release of the Iras data is called Iris (Improved 
Reprocessing of the iRAS t/a^«)(Miville-Deschenes, 2005) and has 
been built with a better destroying, a better subtraction of the 
zodiacal light and a calibration and a zero level compatible with 
the far infrared instrument, FIRAS, of CORE. The IRIS maps 
were also smoothed down to a resolution of one degree and 
downgraded to Nside = 64. 

In order to avoid the contamination from the CMB at inter- 
mediate frequencies, 30-200 GHz, we have restricted our study 
to the Galactic plane where the Galactic emissions dominate over 
the cosmological CMB emission. In practice, we selected those 
regions in the Archeops 353 GHz map with intensity above 3000 
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TABLE II 

Range of values considered for the parameters of the thermal 
dust emissivity model. 



TABLE III 

Spectral index of the free-free emission at the Wmap 
frequencies assuming an electronic temperature of te = 8000 k. 



Parameters 


Range 


Step 




[-1.0,4.0] 


0.02 


Td 


[10.0,37.0J 


0.1 



IiKrj or higher. This corresponds to 1391 pixels at Nside — 64 
in the anticenter region. 



III. Diffuse galactic thermal dust emission 

We first study the electromagnetic and spatial properties of 
the thermal dust diffuse Galactic emission. In order to model 
the intensity of the thermal dust emission, we use a simple grey 
body spectrum of the form 



(1) 



where is the spectral index of the thermal dust emission and 
Td is the dust temperature. 

We used the Archeops and Iris 100 /xm maps to characterize 
the dust thermal emission model. We fitted the data to the model 
pixel by pixel using as free parameters /o, Pd and Td, and the 
following likelihood function 



p2 



where and correspond to the data and model at the pixel 
p within the mask and for the observation frequency v (= 143, 
217, 353, 545 and 3000 GHz), is the 1-cr error bar associated 
to Dl, (3d and Td were explored using an uniformly spaced grid 
as defined in Table ID while Iq was computed via a linear fit 
for each pair (/Sd ,Td). The instrumental noise in the Archeops 
maps has been estimated using simulations of the noise in 
the Archeops Time Ordered Information (TOI) following the 
method described in Macias-Perez et al. (2007). The variance 
per pixel was calculated from 250 simulated noise maps at one 
degree resolution and Nside = 64. The erro r bars for the Iris 
data at 100 /xm were set to 13.5 % following (iMiville-Deschene^. 
1200^ as they are dominated by calibration errors. 



central frequency (GHz) 


23 


33 


41 


61 


94 


/5ff 


-2.090 


-2.093 


-2.095 


-2.099 


-2.103 



TABLE IV 

Range of the parameters considered for the anomalous and 
free-free emission models. 



Parameters 


Ranges 


Step 




[-3.7,-2.3] 


0.01 


Te(K) 


[4000.0, 14000] 


1000 



scope of this study and does not have any consequence in the 
following study. Table U presents the rms of the Archeops and 
Iras data as well as the rms of the residuals after subtraction 
of the dust model. The last column of the table represents the 
mean standard deviation of the noise in the original maps. We 
observe that except for the 5000 GHz data the residuals are of 
the order of magnitude of the noise. 

spectral index anomalous 




■2,3 mKRj 



Fig. 4. Map of the spectral index of the anomalous emission, (right). 



In Figure [T] we present maps of the dust temperature and 

spectral index within the considered mask. We also show the 1^. DiEEUSE GALACTIC EREE-EREE AND SYNCHROTRON 

statistical uncertainties on these parameters. As expected the EMISSIONS 

errors increase significantly on the edges of the maps. These order to estimate the contribution of the diffuse galactic 

noisy pixels will be excluded from the analysis hereafter, free-free emission, which is expected to be important at the 

We can also notice that in the inner regions the statistical Wmap hands, we use the Extinction-Corrected Halpha Foreground 

errors are significantly smaller than the observed dispersion Template (Ha) map built by Finkbeiner (2003). This map was 

for the two parameters. We observe that the mean dust computed using daisL fromihe Virginia Tech Spectral line Survey 

temperature is 20.0 K with 2.1 K dispersion, while the mean (VTSS), for the North and of data from the Southern H-Alpha 

instrumental uncertainties are of the order of 1 K. In the same sky Survey Atlas (SHASSA) for the South sky. Correction factors 

way, the mean dust spectral index is 1.40 with a dispersion of ^re applied to take into account dust absorption iFinkbeiT^ 

0.25, and the mean instrumental uncertainties of the order of 0.1. 2003). We started from a map at resolution Ns^de = 512 and 

downgraded it, as the other maps, at a resolution of Nside — 64. 

Figures ??and|2]compare the best-fit thermal dust model to the order to obtain a template of the free-free emi ssion at 23 

IRAS and ARCHEOPS data. From left to right we show the data, gHz using the Ha map, we follow iBennett et alJ (l2503h . In 

the model and residuals for all frequencies. For the ARCHEOPS antenna temperature units and defining the emission measure 

data the residuals are at most 10 % of the total intensity. In ^s EM — f nldn, one can write 
the case of the Iras data the model reproduce rather well the 

100 /xm map. However, at 60 /xm, the residuals are important //. i^) - ^ aapm [1 + Q.22/n(re/8QQQi^) - 0.14/n(z//41GHz)] 

and the model is not able reproduce the structure in the data. a J - cm-e (zy/41GHz)2(Te/8000i^)i/2 

Residuals can be as important as 60 % of the total intensity. This (3) 

can be expla ined by the presenc e of a hotter dust component as The intensity of the Ha I{R) emission (in Rayleigh units) is 

discussed in iDesert et alJ (Il990l) . This component is out of the defined by 
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SOOOii^ 



-1/2 



0.34/n 



. SOOOii^ / 
(4) 

Thus the intensity of the free-free emission ( in mK rj) is given 
as a function of the intensity of the H-a emission (in Rayleighs) 
by 



T - ^-^^rrm (l + Q-22/n(3^)-Q.14/n( 
^ff- 0.44' ^"^^ 3 



^)) 



(5) 



We have extrapolated this free-free emission template at each 
of the Wmap frequencies assuming that the electromagnetic 
spectrum of the free-fr ee emission is well re presented by a power 
law of the form ui^ff ^Bennett et ^ 120031) 



10.48 + 1.51n(Te/8000i4:) - ln( 



(6) 



41GHz 



We set a sta ndard value for t he electronic temperature at 8000 
K, following (Otte et all 120021) . The values of the spectral index 
obtained at the Wma p frequencies assuming these hypothesis 
are given in Table [Till 

In order to model the synchrotron contribution we used 
the 408 MHz all-sky continuum survey as a template map. We 
extrapolated it at all the considered frequencies assuming a 
power law like electromagnetic spectrum in antenna temperature 
with fix spectral index that we set to -2.7 OBennett et aU y2003h . 



In the second column of Table ?? we present the rms of 
the residuals after subtraction of the Galactic thermal dust, 
synchrotron and free-free emission models. These residuals are 
significant: up to 90 % of the original emission (first column 
of the table). We have observed both point like and diffuse 
structures in these residuals. The former are more probably 
related to uncertainties in the modeling of the free-free emission. 
By contrast, the extra diffuse emission is most probably related 
to anomalous emission. This hypothesis is considered in the 
following section. 

V. Study of the anomalous emission 

In the previous section we concluded that the observed emis- 
sion in the range from 23 to 94 GHz can not be explained only 
by the combination of the canonical Galactic diffuse emission: 
thermal dust, soft synchrotron and free-free. Indeed, we have ob- 
served that in some compact regions there seems to be extra free- 
free emission with respect to the predictions from the Hq, template 
. Furthermore, the diffuse emission is underestimated in general 
indicating either an extra component or a softer synchrotron 
component. In order to investigate these two problems we have 
considered a two component model composed of free-free and 
anomalous emissions in addition to diffuse thermal dust emission. 
We assume that the free-free and the anomalous emissions follow 
a simple power-law model such that 



(7) 



where M^^ are the observed maps in Krj units at the frequency v 
after subtraction of the contribution from thermal dust. Finally, 
we consider 4 free parameters in the model: the normaliza- 
tion coefficients A sync and A//, the spectral index 15 s of the 
anomalous component and the free electron temperature (??). To 
simplify the fitting procedures we vary 15s and Te in the ranges 
shown in Table IIVI Notice that we have not explicitly consider 
the canonical synchrotron emission in this model. Indeed, our so 



called anomalous component will be a mixture of real anomalous 
emission and canonical synchrotron emission. 

We fit this two-component model to the dust subtracted 
Wmap maps and to the 408 MHz map for which the thermal 
dust emission is negligible. As discussed before, the uncertainties 
on the Wmap data have been calculated assuming anisotropic 
white noise on the maps. We compute the variance r pixel using 
the variance per single observation provided on the LAMBDA 
website and maps of the number of hit counts. For the 408 MHz 
map we assume 10 % uncertainties as discussed in Section [III 
It is important to notice that an alternative three component 
component model (including free-free, canonical synchrotron 
and anomalous emission) would imply at least 6 free parameters 
to be fitted on only 6 sky maps. That is why we have chosen to 
consider a two-component model only. 

From the results of the fit we observe that the anomalous 
emission seems to dominate the diffuse component at 1 GHz 
while the free-free emission seems to be mainly located in 
few compact regions. In Figure IH we present the map of the 
reconstructed spectral index for the anomalous emission, /3s. 
We observe that the anomalous emission seems to be well 
represented by a power-law with average spect ra ind e x of - 
2.5. Similar results have been found by Bennet t et alj d2003l) : 
Hinshaw et al. (2007) who claim conclusive evidence for hard 
synchrotron emission. In our analysis we did not dispose of 
data in the frequency range from 10 to 20 GHz to discriminate 
between this hypot hesis and spinning dust emission (refer to 
iDraine & LazarianI |[l998a) for a more complete review on 
spinning dust emission). It is important to notice that currently 
spinning dust emission has mai nly being found in par t icular 
Galactic clouds (see for example | de Oliveira- Cost a et alj dl999L 
'2002!); "Lagache (2003); Hildebrandt et al. (2007); Watson et al] 
(2005); Ysard et al. (2010); Dickinson et al. (2010); Bot et alJ 
(2010); Planc k Collaboration et a l. (2011c) ). Regarding the elec- 
tron temperature, we have found that a physically accessible 
temperature is associated to only 40 out of 1039 pixels considered. 
These pixels corresponds to the intense regions on the free-free 
map at 1 GHz. For the othe r pixels, th e temperature is higher 
than the upper limit allowed dOtte et all 120021) and then can not 
be linked to the free-free emission. 

VI. Summary and conclusions 

We have presented in this paper a detailed analysis of the 
Galactic diffuse emissions at the Galactic anti center in the 
frequency range from 23 to 545 GHz. We have shown that a 
simple grey-body model can be used to describe the thermal 
dust emission in the frequency range from 100 to 3000 GHz. We 
find a mean temperature of 20 K with an intrinsic dispersion of 
2.1 K and a spectral index of 1.4 with intrinsic dispersion of 0.25. 
These values are significantly larger and lower than expected 
from canonical models of the dust emission, T dust ^ 17 K 
and I5dust — 18 — 2 (see f o r exam ple IITnkbeine r et alj dl999l) ; 
iPlanck Collaboration et al J d2011al)). The sa me kind of results 
have been found by iPlanck Co llaborati on et all d2011bl) . 
although as they fixed the spectral index to 1.8 they obtain 
a lower temperature of 14 K. We have performed a similar 
analysis fixing pdust = 18 and we have also obtained lower 
dust temperatures. At high frequencies (above 3000 GHz) extra 
hot thermal dust emission fr om small dust grai ns is needed to 
account for the observations dDesert et all 119901) . 

The former dust model have been used to extrapolate the 
thermal dust emission to microwave frequencies from 23 to 
100 GHz. After subtraction of the thermal dust emission we have 
shown that the microwave data can not be simply explained by a 
combination of free-free and canonical synchrotron emission. A 
more detailed analysis including AME has shown that the latter 
can be well approximated by a power-law of average spectral 
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TABLE V 

RMS OF THE WMAP DATA AND OF THE RESIDUALS AFTER SUBTRACTION OF THE DUST, FREE-FREE AND STANDARD SYNCHROTRON MODEL AND OF THE 
DUST, FREE-FREE AND ANOMALOUS EMISSION MODEL COMPARED TO THE STANDARD DEVIATION OF THE NOISE. 



Frequency (GHz) 


Data mis (mK/^j) 


Residual DFS mis (mKjij) 


Residual DFA mis (mK/^j) 


Noise standard deviation (mK jij) 


23 


2.03831 


2.02387 


0.353262 


0.183557 


33 


0.907433 


0.884012 


0.0881647 


0.0612279 


41 


0.562764 


0.531677 


0.0372230 


0.0438113 


61 


0.256907 


0.205086 


0.0377750 


0.0205909 


94 


0.209367 


0.112932 


0.0530679 


0.0239733 



index —2.5 in Krj units. This anomalous emission seems to 
dominate the diffuse emission at microwave frequencies while 
free-free emission seems to be located in few compact regions. 
Indeed, we have found that outside those regions the data 
required electron temperature has not physically meaningful. 

The spectral index found for the an omalous emission is 
consistent w ith hard synchrotron emission dBennett et all 120031 ; 
lHinshaw~et al., 2007). However, we can not formally con- 
clude on this as our analysis did not include data in the 
1 to 20 GHz that would hel p discriminating thi s hypot hesis 
from spinning dust emission jPraine & LazariaJ (T99Ssl) for 
which conclusive eviden ce hav e been found on some Galacti c 
clouds Ide Oliveira- Cost a et alJ i|1999L I2002D: iLagachd d2003[) : 
Watson et aL (2005) : Ysard et alJd2010D:lDickinson et alJ d 2010D : 



Bot et alJ dioiOh : fPlanck C ollaboration etalJ d2011cD and as 
diffuse emission by [Hildebrandt et alJ d2007D . 
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Fig. 2. Temperature maps (mK^^j) for the Archeops data (left), the thermal dust emission model (center) and residuals (right). From top to bottom we 
present the 545, 353, 217 and 143 GHz maps. 
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Fig. 3. From left to right: Residuals after subtraction of the thermal dust, free-free and anomalous emission models for the 61 (top) and 23 GHz (bottom) 
maps. 



